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<D ! Abstract 

The path integral representation of the transition amplitude for a particle moving in 
curved space has presented unexpected challenges since the introduction of path integrals 
by Feynman fifty years ago. In this paper we discuss and review mode regularization 
of the configuration space path integral, and present a three loop computation of the 
transition amplitude to test with success the consistency of such a regularization. Key 
features of the method are the use of the "Lee- Yang" ghost fields, which guarantee a 
consistent treatment of the non-trivial path integral measure at higher loops, and an 
effective potential specific to mode regularization which arises at the two loop order. We 
also perform the computation of the transition amplitude using the regularization of the 
path integral by time discretization, which also makes use of "Lee- Yang" ghost fields and 
needs its own specific effective potential. This computation is shown to reproduce the 
same final result as the one performed in mode regularization. 
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1 Introduction 



The Schrodinger equation for a particle moving in a curved space with metric (jfu/^x^ 
has many applications ranging from non-relativistic diffusion problems (described by a 
Wick rotated version of the Schrodinger equation) to the relativistic description of particles 
moving in a curved space-time. However it cannot be solved exactly for an arbitrary 
background metric g^x), and one has to resort to some kind of perturbation theory 
A very useful perturbative solution can be obtained by employing a well-known ansatz 
introduced by De Witt |TJ, also known as the heat kernel ansatz. This ansatz makes use of 
a power series expansion in the time of propagation of the particle. The coefficients of the 
power series are then determined iteratively by requiring that the Schrodinger equation 
be satisfied perturbatively. 

Equivalently, the solution of the Schrodinger equation can be represented by a path 
integral, as shown by Feynman fifty years ago 0. One can formally write down the 
path integral for the particle moving in curved space and check that the standard loop 
expansion reproduces the structure of the heat kernel ansatz of De Witt. However the 
proper definition of the path integral in curved space is not straightforward. In fact it 
has presented many challenges due to complications arising from: i) the non-trivial path 
integral measure ||, ii) the proper discretization of the action necessary to regulate the 
path integral. A quite extensive literature has been produced over the years addressing 
especially the latter point 

In this paper we short cut most of the literature and discuss a method of defining 
the path integral by employing mode regularization as it is by now standard in many 
calculations done in quantum field theory. The methods extends the one employed by 
Feynman and Hibbs in discussing mode regularization of the path integral in flat space || . 
It has been introduced and successively refined in || , |7|] and || where quantum mechanics 
was used to compute one loop trace anomalies of certain quantum field theories. The key 
feature is to employ ghost fields to treat the non-trivial path integral measure as part of 
the action, in the spirit of Lee and Yang ||. These ghost fields have been named "Lee- 
Yang" ghosts and allow to take care of the non-trivial path integral measure at higher 
loops in a consistent manner. The path integral is then defined by expanding all fields, 
including the ghosts, in a sine expansion about the classical trajectories and integrating 
over the corresponding Fourier coefficients. The necessary regularization is obtained by 
integrating all Fourier coefficients up to a fixed mode M, which is eventually taken to 
infinity. A drawback of mode regularization is that it doesn't respect general coordinate 
invariance in target space: a particular non-covariant counterterm has to be used in order 
to restore that symmetry ||. General arguments based on power counting (quantum 
mechanics can be thought as a super-renormalizable quantum field theory) plus the fact 
that the correct trace anomalies are obtained by the use of this path integral suggest 
that the mode regularization described above is consistent to any loop order without any 
additional input. 

As usual when dealing with formal constructions, it is a good practice to check with 
explicit calculations the proposed scheme. It is the purpose of this paper to present a 
full three loop computation of the transition amplitude. The result is found to be correct 
since it solves the correct Schrodinger equation at the required loop order. This gives 
a powerful check on the method of mode regularization for quantum mechanical path 
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integrals on curved space. In addition, we present our computation in such a way that 
it can be easily extended and compared to the time discretization method developed in 
refs. ||, which is also based on the use of the Lee- Yang ghosts. This method requires 
its own specific counterterm (also called effective potential) to restore general coordinate 
invariance. As expected both schemes give the same answer. 

The paper is structured as follows. In section 2 we review the method of mode reg- 
ularization and discuss the effective potential specific to this regularization. In section 3 
we present a three loop computation of the transition amplitude. Here we make use of 
general coordinate invariance to select Riemann normal coordinates to simplify an other- 
wise gigantic computation. We check that the result satisfies the Schrodinger equation at 
the correct loop order. In section 4 we extend our computation to the time discretization 
scheme. This is found to compare successfully with the results previously obtained in 
section 3. Finally, in section 5 we present our conclusions and perspectives. In appendix 
A we present a technical section with a list of loop integrals employed in the text. In 
particular, we discuss how to compute them in mode regularization as well as in time 
discretization regularization. 



2 Mode regularization 



The Schrodinger equation for a particle of mass m moving in a .D- dimensional curved 
space with metric g^ u {x) and coupled to a scalar potential V(x) is given by 



where 



H 



dt 



^ + V(x) 



(1) 



(2) 



with V 2 the covariant laplacian acting on scalars. It can be obtained by canonical quan- 
tization of the model described by the classical action 



dt 



m 



g^ u (x)x^x u - V(x) 



(3) 



when ordering ambiguities are fixed by requiring general coordinate invariance in target 
space and requiring in addition that no scalar curvature term be generated by the orderings 
in the quantum potential^. For convenience we will Wick rotate the time variable t — > —it 
and set m = h = 1 to obtain the following heat equation 



9 T 

dt 



l -V 2 + V{x) 



and corresponding euclidean action 

S[x] = 



dt 



-g^{x)x»x v + V(x) 



(4) 



(5) 



1 One could be more general by coupling the particle also to a vector potential (x) . It is simple to 
do so, since mode regularization will respect the corresponding gauge invariance 0. For simplicity we 
set Af_t(x) — in this paper. 



3 



As mentioned in the introduction the heat equation can be solved by the heat kernel 
ansatz of De Witt Q: 



J2 a n (x,y)t n 

n=0 



(6) 



which depends parametrically on the point y^ that specifies the boundary condition 



S u (x-y) 



Here a(x,y) is the so-called Synge world function and corresponds to 



half the squared geodesic distance. The coefficients a n (x,y) are sometimes called Seeley- 
De Witt coefficients^ and are determined by plugging the ansatz (^) into (|]) and matching 
powers of t. 

Now we want to describe in detail how to get the solution of eq. (Q) by the use of a 
path integral which employs the classical action in (|J). Following refs. |6], [7], we write 
the transition amplitude for the particle to propagate from the initial point xf at time ti 



to the final point Xf at time tf as follows 

(x^,tf\xf ,U) = (x^\e l3H \xi) 

where 

S = S[x, a, b, c] = 



x(0)=x f 

T>x exp 



- l -s 

P. 



(It 
i 



g ia/ {x){x^x v + a»a v + lfc v ) + f3 2 (V(x) + V MR (x) 



'MR 

Vx = VxVaVcDc 



up 



(7) 

(8) 

(9) 
(10) 



For commodity we have shifted and rescaled the time parameter in the action, t = tf + (3t 
with (3 = tf — ti, so that — 1 < r < 0. Note that the total time of propagation [3 
plays the role of the Planck constant Ti (which we have already set to one) and counts 
the number of loops. In the loop expansion generated by (3 the potentials V and Vmr 
start contributing only at two loopsQ. The full action S includes terms proportional to 
the Lee- Yang ghosts, namely the commuting ghosts a M and the anticommuting ghosts 
6 M and & 1 . Their effect is to reproduce a formally covariant measure: integrating them 
out produces Vx = n(det g^ u (x(r))) 1 ^ 2 d D x( y r). As we will discuss, mode regularization 
destroys this formal covariance. Nevertheless reparametrization covariance is recovered 
thanks to the effects of the potential Vmr directly included in the action (g). With 
precisely this counterterm the mode regulated path integral in (|7|) solves the equation 



m (jft) in both sets of variables 



:x' 



l -,tf) and 



v Xj , ti 



and with the boundary condition 



{xf 1 , t\xf, t) 

For an arbitrary metric g^ u (x) one is able to calculate the path integral only in a 
perturbative expansion in (3 and in the coordinate displacements £ M about the final point 

2 It is also customary to redefine the a n (x,y) by extracting a common factor Ai(x,y), where A(x,y) 
is a scalar version of the so-called Van Vleck-Morette determinant. 

3 Reintroducing h one can see that the classical potential V must be of order h° while the counterterm 
Vmr is a truly two loop effect of order h 2 . 
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Xj: £ M = xf — Xj. The actual computation starts by parametrizing 

^(r)=<(r)+^(r) (11) 

where ^ (r) is a background trajectory and g M (r) the quantum fluctuations. The back- 
ground trajectory is taken to satisfy the free equations of motion and is a function linear 
in r connecting xf to x^ in the chosen coordinate system, thus enforcing the proper 
boundary conditions 

< g {r) = ^ - CV. (12) 

Note that by free equations of motion we mean the ones arising from (S) by neglecting 
the potentials V + Vmr and by keeping the constant leading term in the expansion of the 
metric ^(x) around the final point a/i, thus making the space flat. Obviously, one could 
also use the exact solution of the classical equations of motion as background trajectory, 
but this would not change the result of the computation. It would correspond just to a 
different parametrization of the space of paths. 

The quantum fields q^{r) in QTTj) should vanish at the time boundaries since the 
boundary conditions are already included in x% At). Therefore they can be expanded in a 
sine series. For the Lee- Yang ghosts we use the same Fourier expansion since the classical 
solutions of their field equations are = = <+ = 0. Hence 

00 

r{r) = £ Csin(Trmr) (13) 

m=l 

where (p stands for all the quantum fields g M , , , <+ . The measure T>x in (110) is now 
properly defined in terms of integration over the Fourier coefficients </>^ as follows 

M D 

Vx = VqVaVbVc = lim A TT TT m dq^da^db^dc 1 ^, (14) 

M— >oo , 

m=l fj.=l 

where A is a constant. Note that this fixes the path integral for a free particle to 

J Vx exp 

where 

Sfree = fjT U^x" + aPc? + IPd*). (16) 

It is well-known that A = (27T/?)~2", however this value can also be deduced later on from 
a consistency requirement. 

The way to implement mode regularization is now quite clear: limiting the integration 
over the number of modes for each field to a finite mode number M gives the natural 
regularization of the path integral. This regularization resolves the ambiguities that show 
up in the continuum limit. 

The perturbative expansion is generated by splitting the action into a quadratic part 
S 2 , which defines the propagators, and an interacting part S int , which gives the vertices. 
We do this splitting by expanding the action about the final point x^ and obtain 

5 = 52 + 5^ = 52 + 53 + 54+... (17) 



i.9 

n 'J free 




A 



£1 

2:i 



(15) 
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win prp 

VV J~LV;±. 






r° 


So — 


/ dr 

f Uj I 




J-i 




,0 


S3 = 


it 










S 4 = 


it 







(18) 
(19) 



d a d p g^{q a q p + - 2g Q £ / V)(^ + (fq v + a'V + \fc v - 2q^ v 



+(3 2 {V + Vj 



MR) 



(20) 



In this expansion all geometrical quantities, like g^ and d^g^, as well as V and Vmr, 
are evaluated at the final point Xp but for notational simplicity we do not exhibit this 
dependence. S 2 is taken as the free part and defines the propagators which are easily 
obtained from the path integral 



(qW(a)) 
<a"(r)a») 
(^(r)c»> 

where A is regulated by the mode cut-off 



-P gTipj) A(r,a) 
P(T(xf) "A(r,tr) 
-2/3 sr{xj) "A(r,a) 



(21) 



A( 



T, <T 



M 

E 

m=l 



7T 2 m 2 



sin 7rmr sin 7rma 



and has the following limiting value for M — > oo 

A(r, (t) = r(a + 1)0(t - a) + a(r + 1)0(<t - r) 



(22) 



(23) 



Note that we indicate *A(r, cr) = J^A(r, cr), A*(t, cr) = ^A(r, a) and so on. Details on 
the properties of these functions are given in appendix A. 
Now, the quantum perturbative expansion reads: 



(xf,tf\xi ,ti) = I Vx exp 



Is 

P. 



(24) 



where the brackets (• ■ ■) denote the averaging with the free action S 2 , and amount to use 
the propagators given in (ffiD in the perturbative expansion. Note that in the last line of 
the above equation we have kept only those terms contributing up to two loops, i.e. up to 
0(f3), by taking into account that £ M ~ 0(/9a), as follows from the exponential appearing 
in the last line of (|4]) after one averages over Note also that having extracted the 
coefficient A together with the exponential of the quadratic action S 2 evaluated on the 
background trajectory implies that the normalization of the left over path integral is such 
that (1) = 1. 
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I 



Using standard Wick contractions and going through a lengthy calculation one gets 



~/ 4 



_Q 2 

2 D 3 



(25) 



1 1 



/312^ 5 V 



2/5' 



dag^dpg 



Xp 



+ P(V + V MR ), 

9o 



(26) 



+Ag ap g^g vX + Ag afl g pA g 

+^U x g up Cf + 2(<fV A - g aX g pP )Ct p + (2g aX g? p - g aP g xp )t p C 



+ (2g^g xp - 
p 9o 



(27) 



This gives the transition amplitude in the two-loop approximation. 

To test its consistency one can use it to evolve in time an arbitrary wave function 

*(x/, tf) = [ d D x i ^g(x i ){x p , t f \x? t U) (28) 



and verify if ty(xf,tf) solves the correct Schrodinger equation^. Since the transition 
amplitude (|24] ) together with the results (p5|) , (|26|) and (|27| ) is given in terms of an 
expansion around the final point (xf,tf), one Taylor expands the wave function \l/(xj,tj) 



and the measure J g(xi) in eq. (|28|) about that point, performs the integration over d 
d ^, and matches the various terms. The leading term fixes A 



d n 

•J- ; 



and the terms of order {3 give 



A =(2*13)-*, 







o. 



(29) 



(30) 



This last equation means that the wave function \l/ satisfies the correct Schrodinger equa- 
tion dD at the final point (x^,tf). 

4 The factor ^/ g{xi) appearing in cq. (]2^) is suggested by the expression of the path integral in 
eq. (Q) which is formally a scalar for Vmr = 0. However general coordinate invariance is broken by 
mode regularization, and recovered thanks to the effects of the counterterm Vmr- Therefore the measure 
appearing in (|2^) should be considered as an ansatz which is verified, for example, by the calculations 
presented in the next chapter. 
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It is interesting to note that the counterterm Vmr appears only in the last line of eq. 
(|26|) . Actually the value of the counterterm reported in eq. (||) has been deduced in || 
by imposing that the transition amplitude would solve eq. (|30|) . General arguments can 
then be used to show that this counterterm should be left unmodified at higher loops. 
In fact one can consider quantum mechanics on curved spaces as a super-renormalizable 
one-dimensional quantum field theory, and check by power counting that all possible 
divergences can only appear at loop order 2 or less in j3. In the next section we are going 
to check that it is so indeed, expelling doubts which have sometimes been raised that 
mode regularization would be inconsistent at higher loops. Thus one can consider mode 
regularization as a viable way of correctly defining the path integral in curved spaces. 

3 The transition amplitude at three loops 

In this section we want to check eq. (|28"D at the next order in (3, which is equivalent 
to showing that the transition amplitude computed by the path integral satisfies the 
Schrodinger equation not only at the point (x^,tf) but in a small neighbourhood of it. 
This computation can be quite lengthy if done in arbitrary coordinates. To make it feasible 
we select a useful set of coordinates: the Riemann normal coordinates centred at the point 
Xj. In such a frame of reference the coordinates of an arbitrary point x M contained in a 
neighbourhood of the origin are given by a vector z^(x) belonging to the tangent space 
at the origin. This vector specifies the unique geodesic connecting the origin to the given 
point x^ in a unit time. In such a frame of reference the coordinates of the origin are 
obviously given by z^(xf) = 0. In what follows we will use Riemann normal coordinates 
which we keep denoting by x^ since no confusion can arise. 

The expansion of the metric around the origin is given by (see e.g. for a derivation) 



Note that the coefficients in this expansion are tensors belonging to the tangent space at 
the origin. This is a property of Riemann normal coordinates. 

In general, the terms contributing to the transition amplitude up to three loops are 
given by 



Clearly the computation would be quite complex in arbitrary coordinates. Fortunately, 
in Riemann normal coordinates many terms are absent since we obtain 




«V(0) + -iV^OKx' 3 + -V 7J R a ^ /3 (0)x a x'V 

+ (^V 7 V 5j R q ^(0) + ^R aiMT pR rv a s (0)\x a ^x s + 0(x 5 ). (31) 




(32) 



S 3 







(33) 
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S 4 
S 6 



= Jjr \^R a ^x a x^(x^x v + a?a v + \fc v ) + (3 2 (V + V MR ) 
= Jjr l^V^Ra^x^x^x^x" + aW + &"c") + (3 2 x a d a (V + V^) 



(34) 
(35) 



dr 
-i 



^VJVsRcc^ + ^R a ^H r , a s)x a x p x' r x 5 (x^ + a^a v + lfc v ) 
4U 45 



+ yX a x /3 a a 9 /3 (\/ + y M i?) 



(36) 



Note that all structures like R. 



V, Vmr and derivatives thereof are evaluated at the 
origin of the Riemann coordinate system, but for notational simplicity we do not indicate 
so explicitly. The computation is still quite lengthy and we get 



2(5 2 4 



s 2 ) = 



\R a pCe h -^Rh- P(V + V MR ), 
o 

-^ y R a pCfC h + ^V a RC h - ^d a (V + V MR )?*, 



(37) 
(38) 



p(i^(y a vpR + vvai?^) + ^R^pR^ee h 

/^VVi^ + ^R^Rf)?? (h ~ h) 

iiRarxRfTW ( J 6 + h) + P 2 V 2 i? + ±Rl + J, 



45 Miy 30 



t Raiiv ^^ 4 J 12 - ^i^fl^W 6 /13 



,4/ "(^ 5| L' 



(39) 
(40) 
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/? 2 



/? 2 



+ T^R^v 2 J14 + -^R 2 alxv p 3 J15, 



72' 



72' 



(41) 



where the integrals I n are listed and evaluated using mode regularization in appendix A. 
Inserting the specific values of the terms arising from the effective potential Vmr when 
evaluated at the origin 



Vmr = -R 

9 a V M R = g^aR 

d a dpVMR = -V a V^i? — —RcxfjLuxRp^ 



(42) 
(43) 
(44) 
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leads us to the following expression for the transition amplitude at the third loop order 



1 



(27T/5) 



T 



1 



2 



12" 



V12 
1 

288' 



24 



360 
1 

"240 



1 



Rapivfi R^ 



1 



1 / 1 



720 



720 

7 1 



R a »R^ + — (— R + V)R, 



12 V12 



la/3 



1 p2 

720 a ^ v P 



480 
1 / 1 



720" 



2 V12 



6 

2 1 



120 



12 



(45) 



This is the complete expression which should be used to test eq. ( |28"D at order (3 2 . A 
straightforward calculation shows that one indeed obtains an identity after making use 
of eq. (|30D. The mode regulated path integral described in the previous section passes 
this consistency check. Therefore it can be considered as a well defined way of computing 
path integrals in curved spaces. 

Before closing this section it may be useful to cast the transition amplitude in a more 
compact form which can be made manifestly symmetric under the exchange of the initial 
and final point. Keeping on using the Riemann normal coordinates (in which we recall 



x. 



and £ M = — x^ = xf) and defining symmetrized quantities as 



A=-[A( Xi ) + A(xf)] 



(46) 



we can write 



-exp 



1 „ 



12 



^(^Ra^R^s + ^V,V s R a p 



360 
1 



R a „uxR^ vX 



720" 
1 



720" 
1 



:Rati,R^ 



^WRo^p + — V a VpR + -V a VpV 



_ p2 

720 aflu P 



1 



1 



1 



-Rl R - — \7 2 R- — V 2 V 
720 aP 120 12 



O(0 



(47) 



From this expression one can extract (by re-expanding part of the exponential) the leading 
terms of the Seeley-De Witt coefficients clq, a%, a,2 for non-coinciding points and obtain, in 
particular, the one loop trace anomalies for the operator H = — |V 2 + V(x) in two and 
four dimensions. 



4 Time discretization 

The computation performed in the previous section was cast in such a way that can be 
easily extended to a different regularization scheme: the time discretization method devel- 
oped in refs. || . Such a regularization was obtained by deriving directly from operatorial 
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methods a discretized version of the path integral. Taking the continuum limit one rec- 
ognizes the action with the proper counterterm, and the rules how to compute Feynman 
graphs. These rules differ in general from the one required by mode regularization. The 
counterterm V w arising in time discretization differs from Vmr, too. 

The time discretization method leads to the following path integral expression of the 
transition amplitude [0] 



(Xj, tf\x^ , ti) 



A 



1 1/4 



g(xi 



where 

Sint 



dr 
i 



\{gA x ) ~ 9^(x f )y^x u + a»a v + Wc v ) + (3 2 {y(x) + V w (x) 



1 



A = (2tt/3)- " 



_D 

' 2 



(48) 

(49) 

(50) 
(51) 



The propagators to be used in the perturbative expansion implied by the brackets on the 
right hand side of eq. Q4"9] ) are the same as in ( [Zip . The only difference is in the prescrip- 
tion how to resolve the ambiguities arising when distributions are multiplied together. 
The prescription imposed by time discretization consists in integrating the Dirac delta 
functions coming form the velocities and the ghosts propagators (thanks to the Lee- Yang 
ghosts they never appear multiplied together) and using consistently the value 0(0) = \ 

for the step function. Note also the presence of the factor [^rj] 1//4 appearing in this 
scheme. 

The result of the calculation has the same structure as the one reported in eqs. ([37]), 
fl38|), (|39p, (^Op, ([HP with the difference that Vmr should be substituted by Vw, leading 
to 



Ir 



v w = 

d a V W = ^aR 



1. 



d a dpV w 



-VcNpR 



1 

24 



(52) 
(53) 
(54) 



and with the following different values of the integrals computed in time discretization 
regularization 



h = 0, J 3 = 0, J 5 = 0, J 10 = 0, /; 



1.3 



1 

12 : 



'15 



1 

12' 



(55) 



The other integrals are as in mode regularization. Inserting all these values back in 
and expanding the coefficient [^rj] 1//4 at the required loop order give the same transition 
amplitude as in (|45|) or, equivalently, in ([47]). Thus this result constitutes a successful 
test on the method developed in 0. 



11 



5 Conclusions 



In this paper we have discussed a proper definition of the configuration space path integral 
for a particle moving in curved spaces. By performing a three loop computation we 
have tested its consistency and checked that one can equally well obtain the perturbative 
solution of the Schrodinger equation by path integrals. This fills a conceptual gap, showing 
that the perturbative description of a quantum particle moving in a curved space obtained 
by De Witt by solving the Schrodinger equation (i.e. using the canonical formulation 
of quantum mechanics [JJ) can equally well be obtained in the path integral approach 
introduced by Feynman fifty years ago. This approach may also have practical applications 
in quantum field theoretical computations when carried out in curved background using 
the world line formalism [ |T0|| . 

We have mainly described the mode regulated path integral. Its definition was ob- 
tained in ||, and [|J by using a pragmatic approach to identify its key elements, and 
needed a strong check to test its foundations. This we have provided in this paper. We 
find that the method of mode regularization is also quite appealing for aesthetic reasons, 
since it is close to the spirit of path integrals that are meant to give a global picture of 
the quantum phenomena. 

On the other hand we have also extended our computation to the time discretiza- 
tion method of defining the path integrals [Q. This method is in some sense closer to 
the local picture given by the differential Schrodinger equation, since one imagines the 
particle propagating by small time steps. It is nevertheless a consistent way of defining 
the path integral, maybe superior at this stage, since one obtains its properties directly 
from canonical methods. As we have seen also this scheme gives the correct result for the 
transition amplitude. 

An annoying property of the two regularization schemes we have been discussing is 
that they both do not respect general coordinate invariance in target space, and require 
specific non-covariant counterterms to restore that symmetry. It would be interesting to 
find a reliable covariant regularization scheme or, at least, a scheme which while breaking 
covariance (e.g. in the decomposition of the action into free and interacting parts) does 
not necessitates non-covariant counterterms. 
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A Appendix 



The function A(r, a) appearing in the propagators is given in mode regularization by 



M 



A(T,<7)= £ 



m=l 



7T 2 n 2 



2 

sin (nrnr) sin (nma) 



(56) 



and leads to the following limiting values as M — > oo, at least in the bulk (we recall that 
-1 < r,a < 0), 

A(t, ct) = r(a + l)0(r - a) + tr(r + 1)6 (a - r) (57) 
•A(r, a) = A A(r, a) = a + 9{r - a) (58) 

A'(r, a) = ^A(r, a) = r + 6 {a - r) (59) 

'A'(t, a) = |:^A(r, <r) = 1 - *(t - a) (60) 

-A(r, a) = A(r, a) = % - r). (61) 

It is also useful to report the following limiting values for coinciding points 

A(r,r) = r 2 + r (62) 
A'(r,r) = 'A(r,r) = r + i. (63) 

Note that at the regulated level one can easily obtain the following identities by inspection 
of (|6|) and its derivatives 

*A*(r, r) + "A (r, r) = <9 r ('A(r, r)) (64) 
A'(t,t) = at r = -1,0 (65) 
9 T (A(r,r)) = 2A'(r,r) (66) 
"A(r,a) = A"(r,a). (67) 

The limiting values given above should be used with care in the perturbative expan- 
sion of the path integral. Rather, one should resort to the proper regularized expressions 
whenever ambiguities arise. In mode regularization we have adopted the following strat- 
egy: one should partially integrate as much as possible to reach expressions which are free 
of ambiguities, and which can be computed directly in the continuum limit. Following 
this procedure we have obtained the following results needed in the text 

(r 2 ('A* + "A) + A - 2r *A)| T = -~ (68) 
(A CA* + "A) -*A 2 )| r = ~ (69) 
(r 3 ("A" + "A) + r A - 2r 2 *A) \ T = \ (70) 
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h = /V(t(-A- + "A)-t-A 2 )| t = 1 (71) 

h = f° dr (r 4 ('A* + "A) + r 2 A - 2r 3 'A) \ T = ~ (72) 

h = J dr (r 2 A ('A* + "A) + A 2 - 2r 'A A) | T = (73) 

h = j\dr(r 2 ACA- + "A)-r 2 -A 2 )| r = -l (74) 

h = jfV(A 2 (*A* + "A)-'A 2 A)| T = ^ (75) 

J 9 = J\r A| T = ~ (76) 

/ 10 = | ^rdcr (V 2 ('A' 2 - "A 2 ) a 2 + 4r 2 *A 2 - 8r 2 'A "A* a 

+ 2 A 2 - 8 A A* a + 4r A *A* cr + 4r *A A* a) = 1 (77) 



r n = y Jdrda (*A* + **A)| T A (*A* + **A)| CT a + r A*| r *A* A*| CT a 

2r ("A" + -A)| T A* A* | CT a + A| T 'A* A| a + A'| T A A* | CT 
— 2A| r *A A"| ff + 2r ('A' + "A)| T A' A| CT + 2r A'| T 'A A'| ff 
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- 2r(-A* + -A)| T AA'| CT -2r A'| r 'A' A\ a aj = -— (78) 

= J jdrda (r 2 ('A' 2 - "A 2 ) A | a + A' 2 A| a - 2r A* 'A* A| a 

+ r 2 'A 2 ('A* + "A)| ff + A 2 ('A* + "A)\ a - 2r A *A ('A* + **A) | a 

— 2r 2 *A 'A' A* \ a - 2 A A* A'| a + 2r A *A* A'| a 

+ 2r 'A A* A"| a ) = l - (79) 
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A 'A A* + 2r A* 'A 2 - 2r A *A "A" ) = — (80) 



'14 



drda ( r A ('A' 2 - "A ) a — r *A A* *A* a + A 'A* 

1 

18 

= J jdrda (A| T ('A' 2 - "A 2 ) A| a - 4 A| T 'A* *A A*| a 
+ 2 A| T *A 2 ('A* + **A)|<j + 2 A*| T A 'A* A*| a + 2 A'| T *A A* A'| a 
- 4 A*| T A *A ('A* + **A) | <j + ('A* + M A)| T A 2 ('A* + "A)| ff ) = ~ (81) 

= J jdrda (A 2 ('A* 2 - "A 2 ) + *A 2 A' 2 - 2 A 'A A* 'A') = ~ (82) 

On the other hand, the time discretization method needs a different prescription in 
order to resolve the ambiguities. It consists in integrating the Dirac delta functions 
whenever they appear (the Lee- Yang ghosts guarantee that they never appear multiplied 
together) and using consistently the value 9(0) = \ for the step function. We present 
now a list of the elementary integrals needed in the text and whose values differ in the 
two regularizations. We have reported both values, the one related to time discretization 
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being included in square brackets. 



jdr t 2 CA- + "A)| T = -i 
Jdr t 3 CA' + -A)| T = i 
jdr T 4 CA- + "A)| T = -A 

//< 
//< 



Idrda r *A A* *A* a = — 



drda A 'A A* 'A* 



r 

.3. 

r 
I. 



r 

.5. 



13" 
"45. 



1 ■ 

18. 



(83) 
(84) 
(85) 
(86) 
(87) 



180 



.360 J 
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